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Abstract. We study a family of correlated one-dimensional random walks 
with a finite memory range M . These walks are extensions of the Taylor's 
walk as investigated by Goldstein, which has a memory range equal to one. 
At each step, with a probability p, the random walker moves either to the 
right or to the left with equal probabilities, or with a probability q — 1 ^ p 
performs a move, which is a stochastic Boolean function of the M previous 
steps. We first derive the most general form of this stochastic Boolean function, 
and study some typical cases which ensure that the average value < i?„ > of the 
walker's location after n steps is zero for all values of n. In each case, using a 
matrix technique, we provide a general method for constructing the generating 
function of the probability distribution of i?„ ; we also establish directly an exact 
analytic expression for the step-step correlations and the variance < i?^ > of 
the walk. From the expression of < i?^ >, which is not straightforward to 
derive from the probability distribution, we show that, for n going to infinity, 
the variance of any of these walks behaves as n, provided p > 0. Moreover, in 
many cases, for a very small fixed value of p, the variance exhibits a crossover 
phenomenon as n increases from a not too large value. The crossover takes 
place for values of n around 1/p. This feature may mimic the existence of a 
non-trivial Hurst exponent, and induce a misleading analysis of numerical data 
issued from mathematical or natural sciences experiments. 
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1 Introduction. 



Soon after its early developments in probability theory [1] the study of random 
walks and their applications has gradually carved its way into the concerns of 
many topics in exact or applied science. In various recent developments, the 
notion of "walk" in the proper sense of the term has gradually blurred, giving 
way to a more general concept where a "step" of the walk is not necessarily 
an algebraic distance between two locations but may become a time interval 
between two consecutive events of a chronological series, or even a sequence of 
real numbers. 

Before mentioning a few significant fields of application, let us quote more 
specifically the notion of "correlated random walk," which underlies the possi- 
bility for a step to depend on the values taken by some (or all) of the previous 
steps. We will address the step correlations and their consequences at length 
in the course of this paper, and will consider therefore the "direct" problem, in 
contradistinction with the usual concern prevailing in applied science which asks 
for answers to the "inverse" problem, i.e., can a given set of data be suitably 
described by a (correlated or uncorrelated) random walk formalism? Examples 
of fields resorting to the inverse problem are neuromuscular [2] and cardiovascu- 
lar medicine [3], structural genetics [4-6], physiology [7], regulation of biological 
rhythms [8]; to these fields pertaining essentially to the medical domain can 
be added the soaring irruption of financial analysis [9-11] in view of market 
predictions. 

In this paper we will focus on the formalism of correlated random walks on 
one-dimensional lattices and discuss the role of memory in some typical models 
which can be solved exactly. With the help of a direct evaluation of the variance 
of the walks, we will be able to prove that all walks with short-range memory 
have a Hurst exponent equal to ^ in the strict limit of an infinite number of steps, 
although the apparent Hurst exponent may look different from i. In particular, 
wc will endeavor to produce a simple argument for a necessary condition yielding 
a Hurst exponent larger than 5. 

2 Correlated walks and memory 

The transition probabilities in a simple random walk on a one-dimensional in- 
finite lattice arc constant and independent of the course of the walk. If these 
transition probabilities depend upon the M previous steps, we have a correlated 
random walk whose memory range is M. Correlated random walks were intro- 
duced by Taylor [12] in an analysis of diffusion by continuous motions. The 
Taylor's model is a one-dimensional random walk on a lattice, in which steps 
can be made to nearest neighbors only. The walker has a probability a to repeat 
his previous move, and a probability 1 — a to move in the opposite direction. 
The resulting correlated random walk has a memory range equal to 1. If a > ^ 
(resp. a < 5) the walk is said to be persistent (resp. antipersistent) . A complete 
analysis of the Taylor's model has been given by Goldstein [13]. More examples 
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of persistent random walks have been studied by many other authors [14-16]. A 
detailed bibliography may be found in Weiss [17]. 

In this paper we study a class of correlated random walks with memory 
range M > 1. Our model may be described as follows. Let (Sk) be an infinite 
sequence of random variables such that, for A: > M, 

{1 with probability 

— 1 with probability ^p 

fM{Sk-i,Sk-2, ■■ ■ , Sk-M) with probability q = 1 - p, 

where /m, called the memory function, is a Boolean function from {—1,1}^ 
onto { — 1,1}. For k < M, the S^s are independent, identically distributed (iid) 
Bernoulli random variables, each Sk taking the values ±1 with equal probabil- 
ities. The values of these M random variables define the initial conditions for 
the walk. The position of the random walker after n steps is 

Rfi = Si + S2 + ■ ■ ■ + Sn- 

The set J^m of all bounded functions defined on {—1,1}*'^ is a Euclidean space 
of dimension 2*^ for the dot product defined by 

^ JM I JM ^ — 
M 

^^fM\^k-l,Sk-2,--- ,Sk-M)fM\Sk-l,Sk-2,--- ,Sk-M), 




where the trace operator tr is a sum over all the possible values of the M 
variables Sk-i, Sk-2, ■ ■ ■ , Sk-M- The set of the following 2*^ functions 

1, 

Sk-l, Sk-2, ■ ■ ■ ,Sk-M, 

Sk-lSk-2i Sk-lSk-3, • • • , Sk-M+lSk-M, 
Sk-lSk-2Sk-3-, • • • ) Sk-M+2Sk-M+lSk-M , 



Sk-lSk-2 ■ ■ ■ Sk- 



M, 



is a complete orthonormal system, and any function in J^m can be written as a 
linear combination of these 2^^ functions. The most general correlated random 
walk with a finite memory of range M is, therefore, defined by a sequence (Sk) 
of random variables such that, for k > M, 



Sk 



1 

-1 

£7, 7 



rSk-jl Sk- 



, Sk- 



with probability 

with probability |p 

with probability ^ajija- •>) 



(1) 



where p + q = 1, 1 < ji < j2 < ■ ■ ■ 
random, are equal either to 1 or to 



< jr < M, and Sjij^—j^, which are non- 
-1. The a's are conditional probabilities 
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(some of them possibly equal to zero) satisfying the completeness relation: 

l<ji<M l<ji<j2<M 

+ X] %u'2-jr H hai2...M = l. 

i<ii<i2<<><M 



As above, for k < M, Sk is a simple Bernoulli random variable which takes 
the values ±1 with equal probabilities. 

In order to simplify this model, we require that, as in the case of a sim- 
ple symmetric random walk, the average value < i?„ > of the position of the 
random walker after n steps should be equal to zero. While this condition is 
automatically satisfied if n < M, a tedious inspection shows that, for n > M, 
the above requirement is satisfied if, and only if, the memory function /m has 
the following simple forms: 

• fM{Sk-i, Sk~2, ■ ■ ■ , Sk-u) is a linear combination of multilinear terms of 

odd degree in the S'fc's; 

• fM{Sk-ii Sk-2, ■ ■ ■ , Sk-Ai) reduces to a single multilinear term containing 
an even number of S'fc's. 

For example, if M = 3, the most general form for the memory function /s, 
satisfying the condition < i?„ >= 0, is either 

/sC^fc-i, S'fc_2, Sfc-s) = aiEiSfc-i + a2e2<S'fc_2 + as^sSfe-a 

+CH23£^123'S'fe_l6'fc_2<S'fe_3, 

or 

faiSk-ii Sk-2, Sk-a) = O'ij^ijSk-iSk-j-, 

with 1 < I < 3 and 1 < j < 3. 

In the sequel we will only consider two cases: either /m is a linear function 
of the M previous steps, or the product of the M previous steps. 

Let us briefly examine these two cases. 

In what follows {Xk) denotes an infinite sequence of iid Bernoulli random 
variables such that, for all positive integers k, we have 

P{Xk = l)=P{Xk = -l) = \, 
1. /m is linear. If {Sk) is an infinite sequence of random variables defined 

by 

(Xk iik<M, 
Sk = \ Xk with probabiHty if fc > M, (2) 
[^EjSk-j with probability quj, if k > M, 
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where, for j = 1, 2, . . . , M, the ej have fixed values chosen a priori equal either 
to 1 or —1, and the Uj are M nonnegative real numbers such that 

ai + 02 H h Am = 1- 

It is clear that, for k < M, < Sk >= 0. Hence 

< Sm+1 > = p< Xm+1 > +gaiei < 5m > H h qaM^M < Si > 

= 

< Sm+2 > = p< Xm+2 > +gaiei < Sm+i > H h qauSM < S2 > 

= 



< SM+e > = p < XM+t > +qai£i < Sm+i-i > H h qaM^M < Se > 

= 0, 

and, therefore, < -R„ >= 0, for all positive integers n. 

2. /m is a single product ofr terms. If (Sk) is an infinite sequence of random 
variables defined by 

{Xk if fc < M, 

Xk with probability pif k > M, (3) 

^hh---3r-Sk-jiSk-j2 ■ ■ ■ Sk-jr with probability q if fc > M, 

where 9 = 1 — p, 1 < ji < i2 < • • • < ir 5J Af, and £jij2---jr is given and equal 
to either 1 or —1. Here again, for k < M, < Sk >= 0. Hence, assuming first 
that r = M, we have 

< Sm+1 > = p < Xm+1 > +9£i2---m < 'S'iS'2 • • • Sm > 

= p< Xm+1 > +qei2-M < Si >< S2 > ■ ■ ■ < Sm > 
= 

< 'Sm+2 > = p < Xm+2 > +q£i2---M < S2S3 ■ ■ ■ SmSm+i > 

= pqei2...M < S2S3 ■ ■ ■ SmXm+1 > 

+q^£i2---M < S1S2SI ■ ■ ■ Sm > 
= q^ei2...M <Si>= 

where in the expression of <S'm+2 we have replaced Sm+i by its expression. More 
generally, 

< Sm+e > = p < XM+e > +q < S^Se+i ■ ■ ■ Sm+c-i > 

= pq < SgSi+i ■ ■ ■ Xm+i-1 > +q^ < S(,-\SjS}_^i ■ ■ ■ S\[_^_^_2 > 
= q^ < Si-x >, 

which shows that for all positive integers k, < Sk >= 0, and, consequently, 

<Rn>=0. 
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Along similar lines, it can be shown that the result < i?„ >= remains valid 
for r < M , although its proof requires a larger number of substitutions. 

In the following sections we will show that, for any finite value of M, we 
can determine exactly the probability distribution of < i?„ >, using a transfer 
matrix method. It will turn out that the rank of the matrix is equal to 2^. Due 
to the increasing complexity of the eigensystem problem with increasing M, 
this method becomes rapidly cumbersome, although it rather easily provides a 
numerical solution if the number of steps n is not too large. In order to illustrate 
the method, we will first briefly review the Goldstein's solution of the Taylor's 
model. 

Since the asymptotic behavior of the variance < i?^ > is of particular in- 
terest to natural scientists, we will directly derive its exact expression through 
a recursive method that applies at any finite space dimensionality. From this 
expression we will be able to show that the asymptotic behavior < i?^ >, when 
n goes to infinity, exhibits a crossover phenomenon. 

3 Taylor's model. 

The Taylor's model corresponds to Model ^ for M — 1, that is, 
and, for all /c > 1, 




Xk , with probability p 

eSk~i, with probability q ~ I — p, 



where (Xk) is a sequence of iid Bernoulli variables such that, for all positive 
integers k, 

P{Xk = -1) = P{Xk = 1) = ^, 

and £ is a non-random number equal to ±1. 

If Pn{k) denotes the probability that the walker will be at site k after the 
nth step has been completed, we have 

Pn{k) =pt{k) +Pn{k), 

where the superscript ± refers to the sign of the n-th step, that is 

p^{k.n) = P((i?„ = fc) n (5„ - ±1)). 

Taking into account the conditional probabilities: 

P(5„ - 1 I ^„_i = 1) = P{Sn = -I I Sn-1 = -1) = \{P+{1+ e)q) 
P{Sn = 1 I Sn-1 = -1) = P{Sn = -1 I ^„-l = 1) = 1{P+{1- e)q) , 
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we obtain the following recursion relations 

P^ik) = l{p+{l + e)q)p+_,{k-l) + ^{p+{l-s)q)p-^,{k-l) 

= \{^ + eq)pt-i{k-l) + l{l~eq)p-_,{k-l) (4) 
Pnik) = lip+il-e)q)p+_,{k + l) + ^ip+{l + e)q)p-_,ik + l) 

= ^{l-sq)p+_,{k + l) + ^{l+sq)p-_,{k+l), (5) 

or, in a more condensed form which will be useful when we consider the case 
M> 1, 

PnW = ^ (1 - eq)p+_,{k - a) + i (1 + eq)p-_,{k - a), (6) 

with CT = ±1. The probability a in Goldstein's notations coincides here with 
i (! + £<?)• 

The generating function of the probability distribution of the random walk 
is defined by 

Mx) = f+ix)+f-{x), 

with 



k 



k 

where the summation index k runs from — n to n with step 2. Let 



in{x) 



and 



Then 



M(a;) 



\{l + £q) X 5 (1 - eq) x 
\{l-eq) i \{l + eq) i 



f„(x)=M(a;)f„_i(a;). (7) 
Iterating this recursion relation, we obtain 

in{x)=yr-\x)h{x), (8) 

where 



^^^^ = f-( \ 
\ fi (x) 
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2x 



Let R(a;) be the transformation matrix such that 

/ Ai(a;) \ 

D(x) = R-^(x) M(x) R(x) = 

where Ai(a;) and X2{x) are the eigenvalues of M.{x). Then 

f„(x) = nix) D"-i(a;) R-^(a;) fi(a;) 

The above relation involves the eigenvalues and eigenvectors of M(.x) which 
can be readily obtained. However, the exact expression for the generating func- 
tion of Pn{k) is cumbersome to manipulate, and does not easily yield the vari- 
ance. Therefore, we resort to an alternative and more direct method, which 
allows for a recursive evaluation of < i?^ >: 

<Rl> = < (5i + 52 + • • • + 5„)2 > 

+2 < 5„ (Si + 52 + • • • + Sn-i) > + <Sl> 
= < Rl_^ > +2 < 5„ (5i + 52 + • • • + Sn-i) > + <Sl>. 



As will appear below, the correlation coefficient < S„S„-k > does not depend 
explicitly upon n, thus we use the notation 

Cfc =< SnSn—k >) 

and obtain 

< Rl >=< Rl_i > +2(ci + C2 + • • • + c„_i) + 1, 
where we have taken into account that < >= 1. Evaluating the c^'s, we get 

Ci = < SnSn-l >=p < XnSn-l > + EQ < 5^_i >= Sq 

C2 = < SnSn-2 >= P < XnSn-2 > + Sq < 5„_i5„_2 >= SqCi = 

Cfe = < SnSn-k >=P < XnSn-k > +£9 < Sn-lSn-k >= ^Ck-l = (e?)''- 

Substituting the c^'s in the expression of < i?^ > we obtain 

<Rl> = <Rl_^> + 2{eq + q^ + --- + {eqr-^) + l 

- <^n-l> + 2 +1, 

and iteration of this recursion relation yields 

<^"> = ^+ (i_^^g)2 {{n-l)eq-nq^ + {eqr+') 

2 

= n 



{l-e{l-p)Y 
((n - 1)£(1 -p)-n{l-pf + (1 - . (9) 
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It is quite remarkable that for a similar random walk on a d-dimensional 
simple cubic lattice, the above result is preserved. Of course, in this case, the 
sequence of one-dimensional X^'s is replaced by a sequence of d-dimensional 
random vectors X^'s, and each elementary step can be performed along each 
of the d axes, either in the positive or the negative direction. For any fixed 
nonzero value of p, when n goes to infinity, the term (1 — goes to zero, 

and < > behaves as 

1+eq 
-j n, 

I — eq 

and the random walk is eventually Gaussian at very large n. However, when q 
goes to 1, according to whether e is equal to +1 or —1, the prefactor of n goes, 
respectively, either to infinity or to zero. Consequently, if p is very small, the 
asymptotic behavior of < Rf^ > should be carefully analyzed, distinguishing the 
two cases e = 1 and e = — 1. 

3.0.1 Persistent random walk (e = 1). 

Replacing e by +1 in |^, and expanding the resulting expression in powers of p, 
we obtain 

<Rn>=n+^[^- + 0{p) 

which shows that, as expected, when p tends to zero, < R^ > tends to n^. 
Therefore, for a small fixed value of p, when n increases, we expect, for the 
variance, a crossover from a behavior to an n behavior. This crossover is 
controlled by the term (1 — in ^, and the above expansion of < R^ > 

shows that the relevant parameter is the variable a = np. To characterize this 
behavior, let define the exponent 

a log < Rl > 

E{n,p) = — , 

alogn 

which is twice the Hurst exponent. A simple calculation yields 

, „ , a(l 



,p \ = 

p^o \p J a + e " — 1 

This expression shows that, if n is large but small compared to l/p (i.e., a is 
small), the exponent E tends to 2, whereas if n is large compared to l/p (a is 
large), E tends to 1. This result is illustrated in Figure |]for p = 0.00005. Note 
that the crossover takes place around n = e^^ ~ 6 10"*, np being of the order of 
unity. 

The crossover from E{n,p) = 2 to E{n,p) = 1, clearly observed in the above 
figure, where p has a small fixed value and n varies, can also be evidenced from 
the figures representing the probability distribution for a fixed value of p: 

• Figure |^: at small p {np ^1) most of the weight of the distribution is 
concentrated at the edges of the distribution range. 
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Figure 1: Exponent E{n,p) of the persistent walk as a function of log n, for 
p = 5 IQ-^ and n,„ax = 1-6 lO'^. 



• Figure ^: at large p {np 3> 1) the significant weight of the distribution 
spreads out in a Gaussian-hke profile, with a width roughly proportional 
to y/n. 



0.Q25 
.*0 . 02 
." 0.015 
0.01 



Figure 2: p„(fc) for M — 1, e — 1, 
p = 0.002 and n = 100. 



Figure 3: p„(fc) for M — 1, e = 1, 
p = 0.2 and n = 100. 



The fulfillment of both conditions (with perhaps the possibility that, in the 
case corresponding to np ^ 1, the localization takes place at various points of 
the distribution range, and not necessarily at the edges, as, for example, for 
the antipersistent walk discussed below) will be from now on considered as a 
sufficient condition for the occurrence of a crossover. 



3.0.2 Antipersistent random viralk (e = — !)• 

Clearly, if p — the walk will be a zigzag walk, and < R^^ >= (resp. < R^^ >= 
1) if n is even (resp. odd). If p is small, a calculation similar to the previous 
one shows that: 
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• If ^ 1, < i?^ > behaves as ^np, and the oscillatory terms related to 
the parity of n are negligible. 

• If np <C 1, we find that < > behaves as np when n is even, and as 1 
when n is odd. Both expressions are the order of unity or less, so that < R^^ > 
behaves as ri°. 

Figure^ represents the exponent characterizing the behavior of the variance. 
In order to avoid, for values of n such that np <C 1, the spurious behavior of 
< R^ > in n for even values of n, we have only considered the odd values of n. 

1 

0.8 
. 6 
0.4 
0.2 



Figure 4: Exponent E{n,p) of the antipersistent walk as a function oflogn (n 
odd), for p — 5 10~^ and rimax — 1-6 10''. 



4 Linear memory function. 

We first address explicitly the case M = 2, as a significant illustration of the 
complexity which manifests as soon as M > 1. Even in this closed form 

for the probability distribution is already difficult to write down. However, 
we provide a formalism that can be fruitfully applied, at least for numerical 
purposes, to larger values of M. 

4.1 Probability distribution for M = 2. 

Let p'^i°'^(fc,n) denote the probability that, after n steps, the random walker is 
at site k and the steps Sn-i and 5„ are , respectively, equal to cri and (T2, that 
is 

pl^'Hk) = PUR,, = k)n ((5„_i = ai) n (5„ = aa))). 
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If Sn~i — cTi and 5„ — (J2, then i?„_i = fc — (T2, and p^i'^^(fc) is a linear 
combination of two probabilities of the form p,^'^\(fc — (72) where cr = ±1: 

= \{^ + q{ai(yia2ei+a2(J2e2))pll\{k-a2) 

+ ^ (1 + q{ai<yi02£i - a20'2£2)) P,^-U^ ^ ^2)- (10) 

To this equation we associate the equation expressing p'^^"^ (k) , where 0^2 = —0-2, 
since this probability is a linear combination of the two probabilities — 
02) = p'^'^i{k + (T2) where cr = 1 or —1. We have 

Pn'^^'ik) = lil + q{ai<Jia2ei + a2a2e2))p+Z\{k-a2) 

+ ^ (1 + 9(aicriCT2ei ~ a2CT2e2))p,T-i(fc ~ ^2)- (11) 
The generating function of the probability distribution of < R„ > is defined by 

fn{x) = /++(X) + (X) + f-+ix) + J--{xl 

with 

rn"''{x)^Y.Pn'''^k) x\ 
k 

where the summation index k runs from — n to n with step 2. From Equations ^ 
and |ll|, replacing ui and tT2 by all their possible respective values, we obtain 



/r(^) 



M^ 



/++i(-^) 



where 

M+(X) : 

and 

where 

M-{x) 



i (1 + g(ai£i + 02^2)) X 5(1 + g(ai£i - a2e2)) x 
i (1 - g(aiei + a2£2)) ^ 5 (1 - g(aiei - 02^2)) ^ 



fn^ix) 

fn^ix) 



M 



.C-iix) 

fuMx) 



i (1 - q{aiei - 0262)) a; i (1 ^ 9(aiei + 0262)) a; 
i (1 + g(ai£i - a2e2)) ^ ^ (1 + g(ai£i + 0222)) ^ 



(12) 



(13) 



(14) 



(15) 



Equations and |IJ can be grouped together in a unique equation that 
reads 

\ / /,t+i(^) \ 



/r(^) 
/„-+(a^) 

V fn-ix) I V 











/„"+i(^) 
/,ti(^) 

V fn-l{x) J 
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Since the column vector on the right hand side of the above equation is not 
correctly ordered, we have to introduce a 4 x 4 permutation matrix, and write 
this equation under the final form 



(16) 



where 



and 



fn(a;) 



M(x) 




/ 1 \ 

10 

10 

\ 1 / 

Iterating Relation |l^, we obtain 



f„(a:) = (M{x) Vi 



h{x), 



(17) 



where 



f2(x) = 



V f2-{^) J 



I 4 \ 



Further progress toward the expression of the generating function /„ (x) re- 
quires the explicit diagonalization of the 4x4 matrix M(x)'p4, followed by 
the procedure outlined in the M = 1 case. The formal derivation of fn{x) 
is extremely complicated, and would reveal untractable for M > 2 as will be 
seen below. However, the formalism is well suited to numerically evaluate the 
probability distribution of i?„ recursively using |l^. Figures ^ and |6| show the 
probability distribution pn{k) when ai — a2 = 0.5 and £i = £2 = 1 for two 
values of p to illustrate cases <C 1 and np^ 1. 

Clearly, at the light of the comment at the end of subsection 3.0.1, the set 
of these two figures is a sure sign of the existence of a crossover. 
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Figure 5: p„(fc) for M = 2,a\ = 0.5, 
02 = 0.5, £i = 1, £2 = 1, P = 0.005, 
and n = 50. 



Figure 6: p„(fc) for M = 2, ai = 0.5, 

02 = 0.5, £1 = 1, £2 = 1, P = 0.5, 

and n = 50. 



4.2 Variance for M = 2. 

Correlations are more tedious to evaluate than in the case M = 1. Let 

ci,e =< SeSi+i >; 

then 

ci,i = < S1S2 >= 

ci,2 = < S2SS >=< S2{pXz + q{aieiS2 + a2e2Si)) >= qaiEi 
and, for ^ > 3, 

ci,£ = < >=< Se,{pXe,+i + qaiSiSt + qa2S2St-i) > 

= qaiei + qa2S2Ci^e-i. 

Hence 



Similarly, if 



ci,e = qaisi — ; > 1) 

1 - g'a2£2 

C2,e =< SiSi+2 >; 



(18) 



then, for all positive integers 

C2,t = < SeSe+2 >=< Se{pXi+2 + q{ai£iSe+i + a2£2'S'^)) > 



so that 



qaisici^i + qa2£2 



C2,e = i^^l£lll((l - g02£2)'-^) + (Za2£2 {^ > !)• 
1 — 502 £2 



(19) 



More generally, for r > 2, 

Cr,e = < SeSe+r > = < S(,{jpX(,+r + q{aieiSk+r-l + a2S2Sk+r-2)) > 

= q{a\e\Cr-\,e + a2£2Cr-2,e)- (20) 
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From this recursion relation, we obtain the following expression of c^/ 

^•^■^ " ~J ^ ^1 + ~J ^ ^2 > (21) 

where Ai and A2 are the roots of 

A^ — qaiEiX — ga2e2 — 0. 

It is easy to verify that |Ai| < 1 and IA2I < 1, the equal sign being possible 
if p = 0. For p > 0, these conditions imply that the correlations Cr/ decrease 
exponentially in r, and this, in turn, will ensure a linear dependence in n for 
the variance. 

Since the variance is given by 

<Rl> = <{Si+S2 + --- + Snf > 
n—1 n—£ 

= n + 2 Cr,£, 

1=1 r=l 

to find its its explicit expression, we need to evaluate the double summation. 
Summing first over r, we obtain 

(A2 - Ai)(l - Ai) V 1 ^ qa2e2 



r^l 



(A2 - Ai)(l — A2) V 1 ~ '?a2e2 

where we have replaced C2j by its expression given by Summing now over 
i finally yields 

2^ 2(n-l) f qaiei{X2 - qaiEi) 

<Rn> = n+ — — — — ga2e2 

(A2 - Ai)(l — Ai) \ 1 - ga2e2 

2 gai£i(A2 - gai£i) 1 - {qa2£2Y' 

X ~ X 



(A2 — Ai)(l — Ai) 1 — ga2e2 1 — ga2e2 

2 / gai£i(A2-gai£i) \ Ai(l- A^') 

(A2-Ai)(l-Ai) l-ga2e2 ''"'^7 1 - Ai 

_^ 2 ^ gai£i(A2 - gai£i) ^ Ai(A"~^ - (ga2£2)"~^) 

(A2 — Ai)(l — Ai) 1 — ga2£2 Ai — ga2£2 

2(n - 1) /gai£i(gai£i - Ai) 

' + ga2£2 



(A2 - Ai)(l - A2) V 1 - ga2£2 

2 ^ gai£i(gai£i - Ai) ^ 1 - (ga2£2) 



(A2 — Ai)(l — A2) 1 — ga2£2 1 — ga2£2 

gai£i(gai£i — Ai) \ -'^2(1 — Aj"^) 

gci2£2 ' 



(A2 — Ai)(l — A2) \ 1 — ga2£2 / 1 — A2 
2 ^ gai£i(gai£i - Ai) ^ A2(A2~^ - (ga2£2)"~^) 



(A2 — Ai)(l — A2) 1 — ga2£2 A2 — ga2£2 
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For our purpose, a complete analytic discussion of this expression would 
take us too far. Instead, we will illustrate the behavior of the variance with 
four typical choices for the set of the parameters £i, £2, oi, 02 in Figures ^, ||, ^, 
and |l^. Note that we obtain crossovers, respectively, similar to those obtained 
for the persistent and antipersistent walks of the Taylor's model when £2 = 1 
and £1 is either equal to 1 or —1. On the contrary when £2 = —1, There is no 
crossover, the variance behaves as n for both signs of £1. 





Figure 7: Log-log plot of < > 
vs.n for M = 2, £1 = 1, £2 = 1, 
ai = a2 = 0.5, p = b 10"^, n„iax = 
1.610*. 



Figure 8: Log-log plot of < > 
vs.n for M — 2, ei = 1, £2 — —1, 
ai = a2 = 0.5, p — 
1.610*. 



510-' 





Figure 9: Log-log plot of < R^ > 
vs.n for M — 2, ei = —1, £2 = 1, 



ai 
1.610*. 



02 = 0.5, p = 510- 



Figure 10: Log-log plot of < R^ > 
vs.n for A/ = 2, £1 = — 1, £2 = — 1, 



ai ~ 02 
1.610*. 



0.5, p = 510- 



4.3 Linear memory for M > 2. 

The matrix method used for the case M = 2 is formally easy to generalize. 
It yields, for the generating function, a recursion relation similar to 16, where 
f„ is a 2*^-component vector, M a 2^^ x 2*^ matrix and 7^4 is replaced by 
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a 2 X 2 permutation matrix 'P2M . With increasing vahics of M it becomes 
rapidly untr actable. Unfortunately, due to the large number of parameters when 
M > 2, even the variance < i?^ > is unmanageable. This is not the case for the 
single product memory, as will be seen in the next section. However, there is, at 
least, one case for which we can predict that the variance exhibits the crossover 
observed in the persistent walk of the Taylor's model. This occurs when all the 
Ej {j = 1,2,... , M) arc equal to 1, since in this case, it is easy to verify that 
for p = 0, the variance behaves as n^. 



5 Single product memory function. 

Here again we will first consider the case M = 2 as an illustration of the com- 
plexity which manifests as soon as M > 1. Since we will only consider the case 
of a single product, the notation £jij2 - jV '^i^l be replaced by e. 



5.1 Probability distribution for M — 2. 

As in the linear memory function case, Pn^'^^ (k) denotes the probability P{{Rn = 
k) n {{Sn-i = cri) n {Sn = o'2))), and as above, if Sn-i = cti and S'„ = (J2, then 
Rn-i = k — a2, which implies that p^^'^^(A:) is a linear combination of the 
two probabilities — (72), where a = ±1. Since, with a probability q, 

Sn = eSn-iSn-2, wc havc (72 = saai or a = eaia2- Therefore, using notations 
similar to those of the linear case, we have 





where 



and 



M+{x) 




i(l - qe) X 
1(1 + 9^)^ 



(22) 



where 

i(l - qe) X 
, ^(1 + 9^)^ 

The above equations can be grouped together in a unique equation that 
reads 



M-(a;) = 




(23) 



/ /n++(^) \ 

/n+-(^) 



/M+(x) 
V M-(x) 



/^-i(^) 
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Since the column vector on the right hand side of the above equation is not 
correctly ordered, we have to introduce a 4 x 4 permutation matrix, and write 
this equation under the final form 



(24) 



where 



and 



fn(a;) 



M(x) 



M+{x) 
M-(a;) 



/ 1 \ 

10 

10 

\ 1 / 

Iterating Relation p3, we obtain 



where 



f„{x) = (M{x) ViT^^ f2(x), 



f2(x) = 



/ in-) \ 

V f2-{^) J 



I 4 \ 



V 4x^ I 



(25) 



As in the linear case, further progress toward the expression of the generat- 
ing function jnix) is extremely complicated, and would reveal untractable for 
M > 2. However the formalism is well suited to numerically evaluate the prob- 
ability distribution of i?„ recursively using k3. Figures |l], |l|, |l|, |l| show the 
probability distribution p„ (fc) for two values of p to illustrate cases ^ 1 and 
np ^ 1 when e = ±1. 

At the light of the comment at the end of subsection 3.0.1, Figures |ll|, |l^, 
reveal the existence of crossovers. This property will be confirmed directly 
in the next subsection. 



18 



5.2 Variance for M = 2. 

In order to evaluate recursively the variance, we need to determine the correla- 
tions. We have 



< SeSe+i > = eq < Se{SiSe-i) > 

= eq< Se-i >= 0; 

< SeSi+2 > = eg < Se{Si+iSi) > 

= sq< Se+i >= 0; 

< > = eq < Se{Se+2Se+i) > 



< SeSm > = eq< Se{Sm-iSm-2) > 



= e^q^ < Se{{Sm-2Sm-3)Sm-2 > 



These results provide all the correlations: 



q2im-e)/3 if m - £ is a multiple of 3. 



Note that the correlations do not depend upon e. On the other hand, the 
variance satisfies the recursion relation 



< Rl > = < Rl_^ > +1 + 2 < + S2Sn + ■■■ + Sn-lSn > • 



and, for m > f + 3 




if m — £ is not a multiple of 3 



(26) 
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-40 



-20 



20 
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Figure 11: pn{k) for M = 2, e = 1, 
p = 0.005, and n = 50. 



Figure 12: for M = 2, s = 1 

p = 0.5, and n = 50. 



19 



Figure 13: p„(fc) for M ^2, e ^ -1, Figure 14: p„(fc) for M ^ 2, e ^ -1, 
p = 0.005, and n = 50. p = 0.5, and n — 50. 



Iterating this relation and taking into account |2^, we finally obtain 

<i?2> ^ n + (n - 3)<72 + (n-6)<74 + - •• + («- 3 Ln/3j)g2L»/3J 

^1 _52Ln/3j 



n + nq 



an/3jg2L"/31+2,(^„/3j+l)g2Ln/3j^^ 
(1 - (7^)2 



where [x\ denotes the largest integer less or equal to x. On this expression of 
< i?2 > readily verify that for q = the variance is equal to n, while for 
g = 1 it takes the value 

n+{2n- 3) [n/3j - 3( [n/3j f = \n^/3'\ , 

where \x~\ denotes the smallest integer greater or equal to x. Since the functions 
floor and ceiling are not differentiable, we define the crossover exponent E(n,p) 

by 



E{n,p) = - 



" < Rl+i > (p)- < Rl-i > b) 



2 <Rl>{p) 

Figure |l^ represents the variation of E{n,p) as a function of of logn for p 
510-5. 



5.3 Probability distribution for M > 2. 

The matrix method used for the case M = 2 is formally easy to generalize. 
It yields, for the generating function, a recursion relation similar to 24, where 
f„ is a 2''^-component vector, M a 2*^ x 2*^ matrix and 7^4 is replaced by a 
2M ^ 2^ permutation matrix 'P2M . With increasing values of M it becomes 
rapidly untractable. However, the variance < > can be exactly evaluated as 
we will show in next subsection. 
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10 12 14 16 



Figure 15: Exponent E(n,p) as a function o/logn, for M = 2, p = 5 10 ^ and 
rlrna^ = 1.610'^. 



5.4 Variance for M > 2. 

In order to evaluate recursively the variance, we need to determine the cor- 
relations. By definition, < SiSm >= if £ < m < M. This is also true if 
TO = M + 1: 

< SiSm+1 > = qs < S({SiS2 ■ ■ ■ Sc ■ ■ Sm) > 

= qe < Si X S2 > ■ ■ ■ < Sm >= 0, 

due the fact that the random variables Si,S2,--- ,Sm are independent, and 
< S^ >= 1. ^ 

Considering now the case: < SiSm > (any £ and m > £), wc first note that 

< SeSe+M+i > = qe < Si{S£+iSe+2 ■ ■ ■ S^+m) > 

= {qef <{SiSi+i---St+M-lf > 

= g^ 

whereas, in the general case, 

< SiSe+r > = qe < Se{Si+r-MSi+r-M+l ■ ■ ■ Se+r-l) > 

~ q^ < S^Si+r-iM+l) >, (28) 

where it is assumed that £ + r — {M + 1) is positive but may be smaller than £. 
Iterated application of |2^ in the presence of the initial conditions shows that, 
for m-£^ [f^J + 5 (0 < 5 < M + 1), 

fO if 6 > 

<^,^„>=<^ I (29) 
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Substitution of in the recursion relation satisfied by the variance: 

< Rl > = < Rl_^ > +1 + 2 < SiSn + S2Sn + ■■■ + Sn-lSn > 

yields 

<Rl> = n+{n~iM+l))q^ + {n-2{M + l))q'^ + --- 



n + nq 



M + 1 

1 - q'^lTTTii 



(M + 1)^ g^Lirf^J 



- {M + 



(1 - q'^y 

which is the generalization for any M of the result ^obtained for Af = 2. Note 
that M = 1 is not a particular case of |30| . 

When g = 0, < i?^ >= n, whereas, when p goes to zero {q 1), < R^ > 
tends to 

n 

n+(2n- (M + 1)) 

2 

which is of the order of for large n. These two limiting behaviors are the 
signal of a crossover, which is evidenced in Figure ^ This figure shows that 
the variance is a decreasing function of M at a given n. 
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Figure 16: Log-log plot of the variance as a function of n, for M — 2. 8, 32. 128 
(to he read downwards) , p = 5 10^^ and n„iax = 1-6 10^. 

The general expression of the variance is an even function of q, and, therefore, 
does not depend upon e. Such is not the case for the probability distribution 
Pn{k) of Rn as in Figures ^ and |l|. 

Here the probability distribution are symmetric. This is not always the case 
as can be observed in Figures |ll|, |l2|, |l^, and |l^. 
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Figure 17: p„(fc) for M = 5, e ^ 1, Figure 18: for M = 5, e = -1, 

p = 0, andn = 13. < iifg >= 29. P = 0, andn= 13. < Rj^ >= 29. 



6 Conclusion and perspectives 

We have studied a family of correlated random walks with a finite memory range. 
These walks are a generalization of the Taylor's walk. We have established a 
matrix formalism for the generating functions of the probability distributions 
of these walks. We have also derived in many cases an analytic expression for 
the variance. As a function of the number of steps n, for each of these walks, 
the variance becomes ultimately linear in n in the limit n ^ oo. However, 
at comparatively small n, a different behavior of the variance can be observed 
leading to the existence of a crossover phenomenon. This feature shows that, 
when analyzing experimental results, one finds a Hurst exponent different from 
5, this does not necessarily imply an infinite-range memory. If, experimentally 
feasible, one should increase significantly the number of steps to cross-check the 
stability of the exponent before concluding that the system under consideration 
exhibits a non-Gaussian behavior. 

The ultimate Gaussian behavior of the variance at very large n is a direct 
consequence of the exponential dec;ay of the step-step correlations as a function 
of the absolute value of the difference between step indices. Therefore, in order 
to find a different behavior, one has to consider step-step correlations decreasing 
more slowly than exponentially. This is what may happen in the case of a 
power law decay as can be shown as follows. For a symmetric random walk, the 
variance is given by 

n—l n 

<Rl> = n + 2^ <SeSm> 
e=i m=e+i 

n—l n 

= n 2 ^ ^ < SeSe+r > 

i=l r=l 
n — l n 

= n-f2y^y^cr,£, 

t-l r-l 
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where as above Cr^i =< SiSr+e >■ If we now assume that, for some reason, 

1 

'^^'<^ - ^ 

independently of then 

n-l 



r=l 



„o n — r 



The possible non-Gaussian behavior of < i?^ > can only result from the 
asymptotic behavior at large n of the summation in the right hand side of the 
above relation. According to the value of a, the asymptotic behaviors of 



n — r 



r=l 

depends on a as follows: 

• if a > 1, S{n,a) ^ n(^{a), where C is the Riemann function. 

• if a = 1, S{n, 1) ^ nlogn. 

• if a < 1, S{n,a) - 

Note that the power law decay for the step-step correlations is only a sufficient 
condition to obtain a non-trivial Hurst exponent equal to 1 — a/2. 
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